In [914]:
!pip install arch
In [915]:
import os
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
from datetime import date
from scipy import stats
%matplotlib inline
import warnings
warnings.filterwarnings('ignore')
from datetime import datetime, timedelta
from sklearn.metrics import mean_squared_error as mse
from sklearn.preprocessing import MinMaxScaler
from statsmodels.graphics.tsaplots import plot_acf, plot_pacf
import ipywidgets as widgets
from IPython.display import display
In [2128]:
w = widgets.Dropdown(
options=['SELECT','AAPL', 'ABUS', 'ARDS', 'BABA','BFRI',
'FB', 'GME', 'MCD','PFE', 'PLUG',
'QCOM', 'SENS','TSLA', 'TWTR', 'UUUU'],
value='SELECT',
description ='Stock name:',
)
def on_change(change):
if change['type'] == 'change' and change['name'] == 'value':
print("You have selected %s" % change['new'])
w.observe(on_change)
display(w)
In [2129]:
if(w.value == 'AAPL'):
df = pd.read_csv('/content/Final_AAPL.csv')
if(w.value == 'ABUS'):
df = pd.read_csv('/content/Final_ABUS.csv')
if(w.value == 'ARDS'):
df = pd.read_csv('/content/Final_ARDS.csv')
if(w.value == 'BABA'):
df = pd.read_csv('/content/Final_BABA.csv')
if(w.value == 'BFRI'):
df = pd.read_csv('/content/Final_BFRI.csv')
if(w.value == 'FB'):
df = pd.read_csv('/content/Final_FB.csv')
if(w.value == 'GME'):
df = pd.read_csv('/content/Final_GME.csv')
if(w.value == 'MCD'):
df = pd.read_csv('/content/Final_MCD.csv')
if(w.value == 'PFE'):
df = pd.read_csv('/content/Final_PFE.csv')
if(w.value == 'PLUG'):
df = pd.read_csv('/content/Final_PLUG.csv')
if(w.value == 'QCOM'):
df = pd.read_csv('/content/Final_QCOM.csv')
if(w.value == 'SENS'):
df = pd.read_csv('/content/Final_SENS.csv')
if(w.value == 'TSLA'):
df = pd.read_csv('/content/Final_TSLA.csv')
if(w.value == 'TWTR'):
df = pd.read_csv('/content/Final_TWTR.csv')
if(w.value == 'UUUU'):
df = pd.read_csv('/content/Final_UUUU.csv')
In [2130]:
pd.set_option('display.max_colwidth', None)
In [2131]:
df['Date'] = df['Date'].astype("datetime64[ns]")
In [2132]:
del df['Unnamed: 0']
In [2133]:
df.head(5)
Out[2133]:
In [2134]:
df.info()
In [2135]:
df.shape
Out[2135]:
In [2136]:
sns.set(font_scale=0.8)
In [2137]:
# CHANGE CONTEXT TO poster TO INCREASE FONT SIZES
sns.set_context("talk", font_scale=1.3)
# PLOT OUT BTC-USE'S CLOSING PRICES SINCE 2014
with sns.axes_style("darkgrid"):
fig, ax = plt.subplots(figsize=(18,8))
sns.lineplot(x=df.Date, y=df.Close, color='blue')
ax.set_title('Closing Price')
In [2138]:
# CALCULATE PRICE RETURNS AS DAILY PERCENTAGE CHANGE USING pct_change()
df['returns'] = 100 * df.Close.pct_change().dropna()
In [2139]:
# CALCULATE LOG RETURNS BASED ON ABOVE FORMULA
df['log_returns'] = np.log(df.Close/df.Close.shift(1))
In [2140]:
df.head()
Out[2140]:
In [2141]:
# DROPPING THE 1ST ROW OF DATA
# BECAUSE I SHIFTED IT FORWARD TO CALCULATE RETURNS/LOG RETURNS
df.dropna(inplace=True)
In [2142]:
# PLOT DISTRIBUTION PLOTS OF RETURNS & LOG RETURNS
# AND VISUALLY COMPARE THEM WITH THE STANDARD NORMAL DISTRIBUTION
with sns.axes_style("darkgrid"):
fig, axes = plt.subplots(nrows=2, ncols=2, figsize=(18,12))
axes[0][0].plot(df.returns, color='blue')
axes[0][0].set_title('Returns')
sns.distplot(df.returns, norm_hist=True, fit=stats.norm, color='blue',
bins=50, ax=axes[0][1])
axes[0][1].set_title('Returns')
axes[1][0].plot(df.log_returns, color='green')
axes[1][0].set_title('Log Returns')
sns.distplot(df.log_returns, norm_hist=True, fit=stats.norm, color='green',
bins=50, ax=axes[1][1])
axes[1][1].set_title('Log Returns')
plt.tight_layout()
fig.show();
In [2143]:
# CREATE A FUNCTION THAT CALCULATE REALIZED VOLATILITY
# FROM SAILY LOG RETURNS
def realized_volatility_daily(series_log_return):
"""
Get the daily realized volatility which is calculated as the square root
of sum of squares of log returns within a specific window interval
"""
n = len(series_log_return)
return np.sqrt(np.sum(series_log_return**2)/(n - 1))
In [2144]:
intervals = [7, 30, 60, 180, 365]
vols_df = {}
# ITERATE OVER intervals LIST
for i in intervals:
# GET DAILY LOG RETURNS USING THAT INTERVAL
vols = df.log_returns.rolling(window=i)\
.apply(realized_volatility_daily).values
vols_df[i] = vols
# CONVERT vols_df FROM DICTIONARY TO PANDAS DATAFRAME
vols_df = pd.DataFrame(vols_df, columns=intervals, index=df.index)
In [2145]:
# CHANGING MATPLOTLIB STYLE
plt.style.use(['fivethirtyeight'])
fig, ax = plt.subplots(figsize=(18,7))
for i in intervals:
if i == 7:
alpha = 0.5
lw = 1
else:
alpha = 1.0
lw = 2
ax.plot(vols_df[i], label=f'{i}-Day Interval Realized Volatility',
alpha=alpha, lw=lw)
ax.set_title('Realized Volatility Using Different Interval Windows', fontsize=21)
plt.legend(loc='best', prop={'size': 14})
plt.show();
In [2146]:
INTERVAL_WINDOW = 30
n_future = 7
# GET BACKWARD LOOKING REALIZED VOLATILITY
df['vol_current'] = df.log_returns.rolling(window=INTERVAL_WINDOW)\
.apply(realized_volatility_daily)
# GET FORWARD LOOKING REALIZED VOLATILITY
df['vol_future'] = df.log_returns.shift(-n_future)\
.rolling(window=INTERVAL_WINDOW)\
.apply(realized_volatility_daily)
In [2147]:
df.describe()
Out[2147]:
In [2148]:
df.rename(columns = {'Real_or_Fake_tweet': 'Fake_news'}, inplace = True)
In [2149]:
df = df.fillna(df.median())
In [2150]:
df.isna().sum()
Out[2150]:
In [2151]:
df.info()
In [2152]:
df.shape
Out[2152]:
In [2153]:
df=df.dropna()
In [2154]:
df.dtypes
Out[2154]:
In [2155]:
import matplotlib.pyplot as plt
import seaborn as sns
plt.figure(figsize=(40,15))
sns.heatmap(df.corr(),annot=True)
Out[2155]:
In [2156]:
df.hist(figsize=(20, 32), bins=70, xlabelsize=8, ylabelsize=8);
In [2157]:
df_corr = df.corr()['AvgTrueRange']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with AvgTrueRange:\n{}".format(len(golden_features_list), golden_features_list))
In [2158]:
df_corr = df.corr()['NATR']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with NATR :\n{}".format(len(golden_features_list), golden_features_list))
In [2159]:
df_corr = df.corr()['TRANGE']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with TRANGE:\n{}".format(len(golden_features_list), golden_features_list))
In [2160]:
df_corr = df.corr()['O']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Openness:\n{}".format(len(golden_features_list), golden_features_list))
In [2161]:
df_corr = df.corr()['C']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2162]:
df_corr = df.corr()['E']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2163]:
df_corr = df.corr()['A']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2164]:
df_corr = df.corr()['N']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2165]:
df.columns
Out[2165]:
In [2166]:
df_corr = df.corr()['B5_O_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_O_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2167]:
df_corr = df.corr()['B5_C_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_C_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2168]:
df_corr = df.corr()['B5_E_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_E_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2169]:
df_corr = df.corr()['B5_A_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_A_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2170]:
df_corr = df.corr()['B5_N_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_N_Um:\n{}".format(len(golden_features_list), golden_features_list))
Downward momentum correlation
In [2171]:
df_corr = df.corr()['B5_O_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_O_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2172]:
df_corr = df.corr()['B5_C_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_C_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2173]:
df_corr = df.corr()['B5_E_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_E_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2174]:
df_corr = df.corr()['B5_A_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_A_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2175]:
df_corr = df.corr()['B5_N_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_N_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2176]:
df_corr = df.corr()['Fake_news']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Real_or_Fake_tweet :\n{}".format(len(golden_features_list), golden_features_list))
In [2177]:
df_corr = df.corr()['Downward_momentum_created']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Downward_momentum_created :\n{}".format(len(golden_features_list), golden_features_list))
In [2178]:
df_corr = df.corr()['Upward_momentum_created']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Upward_momentum_created :\n{}".format(len(golden_features_list), golden_features_list))
In [2179]:
df_corr = df.corr()['Verified_status_True']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Verified_status_True :\n{}".format(len(golden_features_list), golden_features_list))
In [2180]:
df_corr = df.corr()['Verified_status_False']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Verified_status_False :\n{}".format(len(golden_features_list), golden_features_list))
In [2181]:
sns.set(font_scale=0.8)
In [2182]:
for i in range(0, len(df.columns), 5):
sns.pairplot(data=df,
x_vars=df.columns[i:i+5],
y_vars=['NATR'])
In [2183]:
df.dtypes
Out[2183]:
In [2184]:
df.isnull().sum()
Out[2184]:
In [2185]:
df.fillna(0, inplace = True)
In [2186]:
df.dropna(inplace=True)
In [2187]:
sns.set(font_scale=0.8)
In [2188]:
corr = df.drop('Close', axis=1).corr()
plt.figure(figsize=(12, 10))
sns.heatmap(corr[(corr >= 0.5) | (corr <= -0.4)],
cmap='YlGnBu', vmax=1.0, vmin=-1.0, linewidths=0.1,
annot=True, annot_kws={"size": 8}, square=True);
In [2189]:
df.describe()
Out[2189]:
In [2190]:
# DROPPING ALL NaN VALUES
df.dropna(inplace=True)
In [2191]:
n_zoom = 365
sns.set_context("talk", font_scale=1.3)
# plt.style.use(['seaborn'])
# VISUALIZE REALIZED CURRENT VS. FUTURE VOLATILITY
with sns.axes_style("whitegrid"):
fig, (ax1, ax2) = plt.subplots(nrows=2, ncols=1, figsize=(18,14))
ax1.plot(df.vol_current, alpha=.8, lw=1, color='gray', ls=':',
label='Current Volatility')
ax1.plot(df.vol_future, lw=1, color='blue',
label=f'Next {n_future} Days Volatility (TARGET)')
ax2.plot(df.vol_current[-n_zoom:], alpha=.8, lw=2, color='gray', ls=':',
label='Current Volatility')
ax2.plot(df.vol_future[-n_zoom:], lw=2, color='blue',
label=f'Next {n_future} Days Volatility (TARGET)')
ax1.title.set_text(f'Future vs. Current Daily Volatility \n Using {INTERVAL_WINDOW}-Day Interval')
ax2.title.set_text(f'Zooming in the Last {n_zoom} Days')
ax1.legend(loc='upper left', prop={'size': 13}, frameon=True)
ax2.legend(loc='upper left', prop={'size': 13}, frameon=True)
plt.tight_layout()
plt.show();
Daily Volatility Distribution¶
In [2192]:
with sns.axes_style("darkgrid"):
fig, ax = plt.subplots(figsize=(10,6))
sns.distplot(df.vol_current, norm_hist=True, fit=stats.norm,
bins=50, ax=ax)
plt.title('Daily Volatility Distribution')
plt.show();
Experiment 2: weekly granularity
In [2194]:
w = widgets.Dropdown(
options=['SELECT','AAPL', 'ABUS', 'ARDS', 'BABA','BFRI',
'FB', 'GME', 'MCD','PFE', 'PLUG',
'QCOM', 'SENS','TSLA', 'TWTR', 'UUUU'],
value='SELECT',
description ='Stock name:',
)
def on_change(change):
if change['type'] == 'change' and change['name'] == 'value':
print("You have selected %s" % change['new'])
w.observe(on_change)
display(w)
In [2195]:
if(w.value == 'AAPL'):
df = pd.read_csv('/content/Final_AAPL.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'ABUS'):
df = pd.read_csv('/content/Final_ABUS.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'ARDS'):
df = pd.read_csv('/content/Final_ARDS.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'BABA'):
df = pd.read_csv('/content/Final_BABA.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'BFRI'):
df = pd.read_csv('/content/Final_BFRI.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'FB'):
df = pd.read_csv('/content/Final_FB.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'GME'):
df = pd.read_csv('/content/Final_GME.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'MCD'):
df = pd.read_csv('/content/Final_MCD.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'PFE'):
df = pd.read_csv('/content/Final_PFE.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'PLUG'):
df = pd.read_csv('/content/Final_PLUG.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'QCOM'):
df = pd.read_csv('/content/Final_QCOM.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'SENS'):
df = pd.read_csv('/content/Final_SENS.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'TSLA'):
df = pd.read_csv('/content/Final_TSLA.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'TWTR'):
df = pd.read_csv('/content/Final_TWTR.csv', parse_dates=['Date'], index_col=['Date'])
if(w.value == 'UUUU'):
df = pd.read_csv('/content/Final_UUUU.csv', parse_dates=['Date'], index_col=['Date'])
In [2196]:
df.columns
Out[2196]:
In [2197]:
df.shape
Out[2197]:
In [2198]:
df.isnull().sum()
Out[2198]:
In [2199]:
df = df.fillna(df.median())
del df['Unnamed: 0']
df.rename(columns = {'Real_or_Fake_tweet': 'Fake_news'}, inplace = True)
In [2200]:
df_weekly = df.resample('W').agg('mean')
In [2201]:
df_weekly.shape
Out[2201]:
In [2202]:
plt.figure(figsize=(40,15))
sns.heatmap(df_weekly.corr(),annot=True)
Out[2202]:
In [2203]:
sns.set(font_scale=0.8)
In [2204]:
df_weekly.hist(figsize=(20, 32), bins=50, xlabelsize=8, ylabelsize=8);
In [2205]:
df_corr = df_weekly.corr()['AvgTrueRange']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with AvgTrueRange:\n{}".format(len(golden_features_list), golden_features_list))
In [2206]:
df_corr = df_weekly.corr()['NATR']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with NATR :\n{}".format(len(golden_features_list), golden_features_list))
In [2207]:
df_corr = df_weekly.corr()['TRANGE']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with TRANGE:\n{}".format(len(golden_features_list), golden_features_list))
In [2208]:
df_corr = df_weekly.corr()['O']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Openness:\n{}".format(len(golden_features_list), golden_features_list))
In [2209]:
df_corr = df_weekly.corr()['C']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2210]:
df_corr = df_weekly.corr()['E']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2211]:
df_corr = df_weekly.corr()['A']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2212]:
df_corr = df_weekly.corr()['N']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with conscientiousness:\n{}".format(len(golden_features_list), golden_features_list))
In [2213]:
df_corr = df_weekly.corr()['B5_O_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_O_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2214]:
df_corr = df_weekly.corr()['B5_C_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_C_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2215]:
df_corr = df_weekly.corr()['B5_E_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_E_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2216]:
df_corr = df_weekly.corr()['B5_A_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_A_Um:\n{}".format(len(golden_features_list), golden_features_list))
In [2217]:
df_corr = df_weekly.corr()['B5_N_Um']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_N_Um:\n{}".format(len(golden_features_list), golden_features_list))
Downward momentum correlation
In [2218]:
df_corr = df_weekly.corr()['B5_O_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_O_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2219]:
df_corr = df_weekly.corr()['B5_C_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_C_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2220]:
df_corr = df_weekly.corr()['B5_E_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_E_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2221]:
df_corr = df_weekly.corr()['B5_A_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_A_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2222]:
df_corr = df_weekly.corr()['B5_N_Dm']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with B5_N_Dm:\n{}".format(len(golden_features_list), golden_features_list))
In [2223]:
df_corr = df_weekly.corr()['Fake_news']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Real_or_Fake_tweet :\n{}".format(len(golden_features_list), golden_features_list))
In [2224]:
df_corr = df_weekly.corr()['Downward_momentum_created']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Downward_momentum_created :\n{}".format(len(golden_features_list), golden_features_list))
In [2225]:
df_corr = df_weekly.corr()['Upward_momentum_created']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Upward_momentum_created :\n{}".format(len(golden_features_list), golden_features_list))
In [2226]:
df_corr = df_weekly.corr()['Verified_status_True']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Verified_status_True :\n{}".format(len(golden_features_list), golden_features_list))
In [2227]:
df_corr = df_weekly.corr()['Verified_status_False']
golden_features_list = df_corr[abs(df_corr) > 0.5].sort_values(ascending=False)
print("There are {} strongly correlated values with Verified_status_False :\n{}".format(len(golden_features_list), golden_features_list))
In [2228]:
sns.set(font_scale=0.8)
In [2229]:
for i in range(0, len(df_weekly.columns), 5):
sns.pairplot(data=df_weekly,
x_vars=df_weekly.columns[i:i+5],
y_vars=['NATR'])
In [2230]:
df_weekly.fillna(0, inplace = True)
In [2231]:
df_weekly.dropna(inplace=True)
In [2232]:
corr = df_weekly.drop('Close', axis=1).corr()
plt.figure(figsize=(12, 10))
sns.heatmap(corr[(corr >= 0.5) | (corr <= -0.4)],
cmap='YlGnBu', vmax=1.0, vmin=-1.0, linewidths=0.1,
annot=True, annot_kws={"size": 8}, square=True);
Weekly volatility distribution
In [2233]:
with sns.axes_style("darkgrid"):
fig, ax = plt.subplots(figsize=(10,6))
sns.distplot(df_weekly.NATR, norm_hist=True, fit=stats.norm,
bins=50, ax=ax)
plt.title('Weekly Volatility Distribution')
plt.show();
